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Abstract. Bona [6J studied the distribution of ascents, plateaux and descents in the class of 
Stirling permutations, introduced by Gessel and Stanley llT4l . Recently, Janson lITSl showed 
the connection between Stirling permutations and plane recursive trees and proved a joint nor- 
mal law for the parameters considered by Bona. Here we will consider generalized Stirling 
permutations extending the earlier results of ||6l, IfTSl . and relate them with certain families 
of generalized plane recursive trees, and also {k + l)-ary increasing trees. We also give two 
different bijections between certain families of increasing trees, which both give as a special 
case a bijection between ternary increasing trees and plane recursive trees. In order to de- 
scribe the (asymptotic) behaviour of the parameters of interests, we study three (generalized) 
Polya urn models using various methods. 



1. Introduction 

Stirling permutations were defined by Gessel and Stanley [[T4l . A Stirling permutation is a 
permutation of the multiset {1, 1, 2, 2, . . . , n, n} such that for each i, 1 < i < n, the elements 
occuring between the two occurences of i are larger than i. The name of these combinatorial 
objects is due to relations with the Stirling numbers, see lfT4l for details. 

Let cr = aia2 ■ ■ •a2n be a Stirling permutation. Let the index i (or the gap + 1)) 
be called an ascent of cr if i = or < Oj+i, let i be called a descent of a if i = 2n 
or > Oj+i, and let i be called a plateau of a if a, = Oj+i. (It is convenient to define 
«o = 0'2n+i = 0; this takes care of the special cases i = and i = 2n.) Note that i runs from 
to 2n, so the total number of ascents, descents and plateaux is 2n + 1. Let Q„ denote the 
set of Stirling permutation of {1, 1, 2, 2, ... , n, n}; we say that these have order n. Bona [6l 
showed that the parameters numbers of ascents, descents and plateaux are equidistributed on 
Qn. Moreover, he showed a central limit theorem for the three parameters. 

A rooted tree of order n with the vertices labelled 1, 2, . . . , n, is an increasing tree if the 
node labelled 1 is distinguished as the root, and for each 2 < k < n, the labels of the nodes 
in the unique path from the root to the node labelled k form an increasing sequence. We 
will consider several families of increasing trees. The first one is the family of increasing 
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plane trees, usually called plane recursive trees, where the children of a node are ordered 
(from left to right, say). Note that plane recursive trees also appear in literature under the 
names plane-oriented recursive trees, heap-ordered trees, and sometimes also as scale-free 
trees. Further families will be defined later. 

Let Tn denote the set of plane recursive trees with n vertices. It was shown by Janson 
[fT8l that plane recursive trees on n + 1 vertices are in bijection with Stirling permutations on 
{1, 1, 2, 2, . . . , n, n}, Tn+i = Qn- Moreover, using this bijective correspondence, he showed 
that the number of descents in the Stirling permutation corresponds to the number of leaves 
in the associated plane recursive tree. Furthermore, using an urn model and general theo- 
rems, see fT5] and also fT6l, Janson showed the joint normality of the parameters ascent, 
descent and plateau. The purpose of this work is to extend this connection between Stirling 
permutations and plane recursive trees in Janson [18], to generalized Stirling permutations. 
In particular, we give a bijection between Stirling permutations on {l'^, 2^^, . . . , ra'^}, where 
here and throughout this work l' := 1, . . . , 1, with I > 1, which we call A;-Stirling permuta- 

i 

tions, and (k + l)-ary increasing trees; moreover we can also relate fc-Stirling permutations 
with a certain family of plane recursive trees, namely A;-plane recursive trees. Concerning 
Stirling permutations of the multiset {l'^, 2^+^, . . . , n'^"'"^}, which we call fc-bundled Stirling 
permutations, we obtain a bijection with certain generalized plane recursive trees, namely 
A; -bundled increasing trees. We also give two different bijections between certain families 
of increasing trees, which both give as a special case a bijection between ternary increasing 
trees and plane oriented increasing trees. Moreover, we will use several different methods, 
combinatorial and probabilistic, to derive several results in this direction. More precisely, in 
order to describe the (asymptotic) behaviour of the parameters of interests, we study three 
(generalized) Poly a urn models. 

The parameter k is fixed throughout the paper, and often omitted from the notation. All 
unspecified limits are as n ^ oo. In the results with a.s. convergence, we assume that the 
random /c-Stirling permutation grows in the natural way by random addition of new labels; 
in the other results, this does not matter. 



2. Preliminaries 

2.1. Generalized Stirling permutations. A straightforward generalization of Stirling per- 
mutations on the multiset {1, 1, 2, 2, . . . , n, n} is to consider permutations of a more general 
multiset {1^^,2^^, . . . ,n^"}, with fcj G N for 1 < i < n. We call a permutation of the 
multiset {1*^^ 2*^^ . ^ j ^ generalized Stirling permutation, if for each i, \ < i < n, the 
elements occurring between two occurrences of i are at least i. (In other words, the elements 
occurring between two consecutive occurrences of i are larger than i.) Such permutations 
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have already previously been considered by Brenti [|71|, |[8l. The number of generalized Stir- 
ling permutations of {I'^i, 2^^^ ^ j^k^j 

n— 1 i 

+ with£, = 5^fc,; (1) 

i=l j=l 

this is easy to see by induction, since the A;„ copies of n have to form a substring, and this 
substring can be inserted in + 1 positions (viz., anywhere, including first or last) in any 
generalized Stirling permutation of {l'^^, 2'^^, . . . ,{n — l)'^"-^}. 

We will consider two cases and give them special names: a k-Stirling permutation of order 
n is a generalized Stirling permutation of the multiset {l'^, 2^^, . . . , ra'^}, and a k-bundled Stir- 
ling permutation is a generalized Stirling permutation of the multiset {l'^, 2^^+^, . . . , 
Here k > 1, but note that 1-Stirling permutations are just ordinary permutations so we will 
usually consider fc-Stirling permutations for k > 2 only; the case k = 2 yields the ordinary 
Stirling permutations defined by Gessel and Stanley [[T4| . 

What we call A;-Stirling permutations was suggested by Gessel and Stanley [fT4l and has 
been studied by Park [i22l l23l l24l under the name A;-multipermutations. 

In the following, let Q„ = Qnik) denote the set of /c-Stirling permutations of order n and 
let Qn = Qnik) denote the number | Qn{k) \ of them. By ([T]), 

Qnik) = \Qn{k)\ = llikz + 1) = r ^^^+^y^ (2) 

For k = 2 this number is just (5n(2) = {2n — 1)!!. In the case A; = 3, we have for example 
one permutation of order 1: 111; four permutations of order 2: 111222, 112221, 122211, 
222111; etc. 

Similarly, let Q„ = Qn{k) denote the set of fc-bundled Stirling permutations of order n 
and let Q„ = Qnik) denote the number of them. We have, by ([U), 

Qn = mk)\ = ii('{k + 2) - 1) = (A: + ^t-'yiBti^t^- 

We define ascents, descents and plateaux of a generalized Stirling permutation o = 
0-10.2 ■ ■ ■ Oi of {l'^^ 2^^^, . . . , n'^"} (where the length i = J2i as before: we let ao = 
a^+i = and say that an index i, with < ? < £, is an ascent, descent or plateau if 

< Oj+i, Oj > Oj+i or ai = Oj+i, respectively. Note that the total number of them is £ + 1. 

We introduce a natural refinement of ascents, descents and plateaux, namely j -ascents, 
j -descents, and j -plateaux. An index i, with 1 < i < i is called a j -ascent, if i is an ascent 
and there are exactly j — I indices i' < i such that Oj/ = a^; i.e., Oj is the jth occurrence of 
the symbol ai, and similarly for plateaux. For a descent i, ai is always the last occurence of 
that symbol (just as for an ascent, a^+i is the first of its kind), and we define a j-descent as a 
descent i < £ such that Oj+i is the the jth occurrence of that symbol. (Note that we choose 
not to allow i = Oori = im these definitions.) 
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Thus, for a generalized Stirling permutation of 2^^^, . . . , n''"}, the possible values of 
j ranges from 1 to maxj ki for j-ascents and j-descents, and from 1 to maxj kt — 1 for j- 
plateaux. In particular, for /c-Stirling permutations, 1 < j < A; for j-ascents and j-descents, 
and I < j < k — 1 for j -plateaux. Note also that if we reflect a /c-Stirling permutation, we get 
a new /c-Stirling permutation, and j-ascents in one of them correspond to (/c + 1 — j)-descents 
in the other. 

Example 1. Consider the 3-Stirling permutation a = 112233321: Index 1 is a 1 -plateau, 
index 2 is a 2-ascent, index 3 is a 1 -plateau, index 4 is a 2-ascent, index 5 is a 1 -plateau, 
index 6 is a 2-plateau, index 7 is a 3-descent, and index 8 is a 3-descent. (Indices and 9 are 
not classified in this way.) 

We are interested in the (joint) distributions of the random variables Xnj, Ynj and Znj, 
defined as the numbers of j-ascents, j-descents and j-plateaux, respectively, in a random 
A;-Stirling permutation (chosen uniformly in Q„(A;)). Note that these trivially are unless 

1 < j < kfor Xn,j and Ynj, and 1 < j < A; — 1 for Zn,j, and that 

k k-1 
^{XnJ + Yn,j) + ^ ZnJ = ku - 1 
j=l j=l 

We further let X„, F„ and Z„ denote the total numbers of ascents, descents and plateaux, 
respectively. Note that, recalling the special definitions at the endpoints, 

k 

Xn = Y,^n,j + l, (4) 

k 

Yn = J2Ynj + l, (5) 
j=i 

k 

Zn = '^ZnJ. (6) 

It is easy to see that a j -ascent with j < k corresponds to a later (j + 1) -descent, and 
conversely, so 

Xnj = Yn,j+U l<J<k-l, (7) 

see also Theorem [21 However, there is no corresponding relation for fc-ascents, of for 1- 
descents, and the total numbers of ascents and descents are typically different, even in the 
case k = 2. Further, since only the last copy of a label can be a descent, 

X„j + Zn,j = n, 1 < J < A; - 1, (8) 

and, similarly or by ([8]), 

+ = n, 2<j<k. (9) 
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Moreover, we are also interested in the distribution of the number of blocks in a random k- 
Stirling permutation of order n. A block in a generalized Stirling permutation o = ai - ■ ■ 
is a substring ap - ■ ■ with = that is maximal, i.e. not contained in any larger such 
substring. There is obviously at most one block for every j = 1, . . . ,n, extending from the 
first occurrence of j to the last; we say that j forms a block when this substring really is 
a block, i.e. when it is not contained in a string i - ■ - i for some i < j. In particular, in a 
/c-Stirling permutation, j forms a block if for any i with 1 < i < j — 1, there do not exist 
indices itlq, . . .rrii^^i , with 1 < niQ < ■ ■ ■ < m^+i < kn, such that cXm^^ = (Tmk+i = 
and cTmj = . . . = a^n^ =3- It is easily seen by induction that any generalized Stirling 
permutation has a unique decomposition as a sequence of its blocks. Note that if we add a 
string [n + to a generalized Stirling permutation, this string will either be swallowed 

by one of the existing blocks, or form a block on its own; the latter happens when it is added 
first, last, or in a gap between two blocks. 

Example 2. The 3-Stirling permutation a = 112233321445554666, has block decomposi- 
tion [112233321] [445554] [666]. 

One may also consider the similar problems for A; -bundled Stirling permutations; simi- 
larly defining random variables Xnj, Ynj and Znj. However, for most results we restrict 
ourselves to /c-Stirling permutations. 

2.2. Generalized plane recursive trees and rf-ary increasing trees. In order to relate the 
fc-Stirling permutations to families of increasing trees we use a general setting based on 
earlier considerations of Bergeron et al. lO and Panholzer and Prodinger [21]. 

For a given degree- weight sequence {(pk)k>o, the corresponding degree- weight generating 
function (f{t) is defined by (^{t) := '^^y^'^kt^ ■ The simple family of increasing trees T 
associated with a degree- weight generating function '^{t), can be described by the formal 
recursive equation 

T = ®x [ipQ ■ {e} u ■ r Uip2-T *r uv'g-r *t *r u ■■■^ = ®x ip{r), (lo) 

where denotes the node labelled by 1, x the cartesian product, U the disjoint union, * 
the partition product for labelled objects, and v^(T) the substituted structure (see e. g., the 
books ||30l , |fT3l ). This means that the elements of T are increasing plane trees, and that a 
tree with (out-)degrees (ii , . . . , (i„ is given weight YYi Vdi ■ By a random tree of order n from 
the family T, we mean a tree of order n chosen randomly with probabilities proportional to 
the weights. 

Let Tn be the total weight of all such trees of order n. It follows from (flOl) that the expo- 
nential generating function T{z) := X]n>i ^"f^f ^he total weights satisfies the autonomous 
first order differential equation 



T\z) = ^{Tiz)), r(o)=o. 



(11) 
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The families that we will consider have degree- weights of one of the two following forms, 
studied by Panholzer and Prodinger 11211 : 



(1+^) 

a 



£r— ;-, for y^o > 0, < — C2 < Ci, generalized plane recursive trees, 



(/?ofl + —1 , for (y9o, C2 > 0, (i := — + 1 G M \ {1}, d-ary increasing trees. 

(12) 

Consequently, by solving (fTTI) . we obtain exponential generating function T{z) 



Vo i 1 



generalized plane recursive trees, 



C2 

and the total weights T„, 



nz) = { ^ . (13) 

^ ' 1 , d-dxy increasing trees. 



(l-(d-l)C22)d-l 



T^ = V,cr\n~l)\{ (14) 

Note that changing Lp^ to ah^ipk for some positive constants a and h will affect the weights of 
all trees of a given order n by the same factor a^lf^^, which does not affect the distribution 
of a random tree from the family. Hence, when considering random trees from these two 
classes, is irrelevant and ci and C2 are relevant only through the ratio Ci/c2- (We may 
thus, if we like, normalize tp^ = \ and either ci or |c2|, but not both.) 

As shown by Panholzer and Prodinger f2Tj, random trees in the two classes of families 
given in (fT2l) can be grown as an evolution process in the following way. The process, 
evolving in discrete time, starts with the root labelled by 1. At step i + 1 the node with label 
2 + 1 is attached to any previous node v (with out-degree d{v)) of the already grown tree of 
order i with probabilities p{v) given by 

d{v)+a Q,;=_i_£i>o, generalized plane recursive trees, 

{d~i)ili ' '^"^y increasing trees. 

Moreover, Panholzer and Prodinger [[2T]| showed that there are only three classes of simple 
families that can be grown in this way (for suitable p(f )): the two classes given in (fT2)) and 
the recursive trees given by Lp(t) = lpqc'^^^^'-^^ with lpq, ci > (which can be regarded as a 
limiting case of any of the two classes above, letting C2 0.) 

Example 3. Plane recursive trees are plane increasing trees such that all node degrees are 
allowed, with all trees having weight 1. Thus (p^ = 1 and the degree- weight generating 
function is ip{t) = j^, which is of the form in (fT2l) with cpQ = 1, ci = 2 and C2 = —1. We 
have 



T{z) = 1 - Vl - 2z, and T„ = 1 ■ 3 ■ 5 ■ ■ ■ (2n - 3) = (2n - 3)!!, for n>l. 
Furthermore, « = — 1 — — = 1, and consequently, the probability attaching to node v at step 
z + 1 is given by p{v) = '^2i-i ■ 
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Example 4. For an integer d > 2, d-ary increasing trees are increasing trees where each node 
has d (labelled) positions for children. Thus, only outdegrees 0, . . . , are allowed; moreover, 
for a node with k children in given order, there is thus (^) ways to attach them. Hence, this 
family is given by vertex weights ipk = (^) and thus the degree-weight generating function 
ip{t) = 1 + which is of the form in with ip^ = I, ci = d - I and C2 = 1. By (fT3l) . 

T{z) = {l-{d-l)zY"^''-''^ -I. 

3. Increasing trees associated to generalized Stirling permutations 

3.1. [k + l)-ary increasing trees, fc-plane recursive trees and A;-Stirling permutations. 

Recall from Example |4] that, for k >1, the degree- weight generating function of (A; + l)-ary 
increasing trees is given by (p{t) = {l + t^^^, i.e. (p^ = 1, Ci = k and C2 = 1. Consequently, 
the generating function T{z) and the numbers T„ of (fc + l)-ary trees of order n are given by 

1 " 
[l-kz)k 

and the the probability of attaching to node v at step z + 1 is given by p{v) = ■ 

Note that T„ = Qn, the number of /c-Stirling permutation, which makes the following 
theorem reasonable. 

Theorem 1 (Gessel). Let k > 1. The family An = An{k + 1) of{k + l)-ary increasing trees 
of order n is in a natural bijection with k-Stirling permutations, An{k + 1) = Qn{k). 

Remark 1. The authors independently derived the result above, and later discovered the 
work of Park [122]. in which Gessel's result was mentioned but the proof only sketched. The 
result of Gessel never appeared in print except this mentioning in Park ^22], to the best of the 
authors' knowledge. We will give a detailed proof of the result above, which has interesting 
consequences regarding the (refined) parameters ascents, descents and plataeux, and also 
number of blocks, which we will state in Theorem [2l 

Remark 2. For /c = 1 we obtain a bijection between 1-Stirling permutations (ordinary per- 
mutations) and binary increasing trees, which is very well known. 

Proof. We use a slightly modified bijection to the one given by Janson in [fTSlI for Stirling 
permutation and plane recursive tree, and use a depth-first walk. The depth-first walk of a 
rooted (plane) tree starts at the root, goes first to the leftmost child of the root, explores that 
branch (recursively, using the same rules), returns to the root, and continues with the next 
child of the root, until there are no more children left. We think of (k + l)-ary increasing 
trees, where the empty places are represented by "exterior nodes". Hence, at any time, any 
(interior) node has k + 1 children, some of which may be exterior nodes. Between these 
k edges going out from a node labelled v, we place k integers v. (Exterior nodes have no 
children and no labels.) Now we perform the depth-first walk and code the (k + l)-ary 
increasing tree by the sequence of the labels visited as we go around the tree (one may think 
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of actually going around the tree like drawing the contour). In other words, we add label v 
to the code the k first times we return to node v, but not the first time we arrive there or the 
last time we return. A (A; + l)-ary increasing tree of order 1 is encoded by l'^. A (A; + l)-ary 
increasing tree of order n is encoded by a string of k ■ n integers, where each of the labels 
1, . . . ,n appears exactly k times. In other words, the code is a permutation of the multiset 
{l'^, 2^^, . . . , n'^}. Note that for each i, 1 < i < n, the elements occurring between the two 
occurrences of i are larger than i, since we can only visit nodes with higher labels. Hence the 
code is a fc-Stirling permutation. Moreover, adding a new node n + 1 at one of the kn + 1 free 
positions (i.e., the positions occupied by exterior nodes) corresponds to inserting the A;-tuple 
(n + I)''' in the code at one of kn + 1 gaps; note (e.g., by induction) that there is a bijection 
between exterior nodes in the tree and gaps in the code. This shows that the code determines 
the {k + l)-ary increasing tree uniquely and that the coding is a bijection. See Figure [T] for 
an illustration. 

The inverse, starting with a fc-Stirling permutation a of order n and constructing the cor- 
responding {k + l)-ary increasing tree can be described as follows. We proceed recursively 
starting at step one by decomposing the permutation as a = aila2l . . . cTfclcTfc+i, where (af- 
ter a proper relabelling) the (jj's are again A;-Stirling permutations. Now the smallest label in 
each (Tj is attached to the root node labelled 1. We recursively apply this procedure to each 
CTj to obtain the tree representation. □ 



Figure 1. The three ternary trees of order 2 encoded by 2211, 1221 and 
1 122; an order 3 ternary increasing tree encoded by the sequence 23321 1. 

Now we relate the distribution of j-ascents, j-descents and j-plateaux in A;-Stirling permu- 
tations with certain parameters in (k + l)-ary increasing trees. In order to do so we introduce 
two kinds of parameters. The parameter Dnj, standing for "jth children" , counts the num- 
ber of nodes in a random (k + l)-ary increasing tree of order n that are the jth children of 
their respective parents, going from left to right, with 1 < j < k + 1. Similarly, the param- 
eter Lnj, Standing for "leaves" of type j, counts the number of exterior nodes that are jth 
children of their parents, 1 < j < k + 1. We thus have 




Ln,j = n- DnJ, I < j < k + 1, 



(15) 
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and, counting the total numbers of interior and exterior children, 

fe+l k+l 

Dn,j =n-l, ^ Ln,j = kn + 1. (16) 

Concerning the number of blocks in A;-Stirling permutation, we introduce one more pa- 
rameter in (k + l)-ary increasing trees. Let LRn denote the number of (interior) nodes that 
have the property that the path to the root consists exclusively of the leftmost or rightmost 
possible edge at each node, i.e., the edge in position 1 or A; + 1, and no other "inner" edges. 
Subsequently, we will call such nodes left-right nodes. The root is trivially a left-right node. 

Theorem 2. Let k > 1. Under the bijection in TheoremU} the numbers of j -ascents Xn.j, 
j -descents Ynj and j -plateaux Znj in a k-Stirling permutation of order n coincide with the 
(shifted) numbers of j -children Dnj, and j -leaves Lnj in a {k + l)-ary increasing tree of 
order n by the formulas 

Xn,j = Dn,j+1, i l£ j ^ k, 

Zn,j = Ln,j+1 = n- Dn,j+1, 1 < j < - 1- 

As a consequence, for the numbers of ascents, descents and plateaux, 

Xn = n — Dn,l = Ln,l, 

= ^ ~ Dn,k+1 = Ln^k+li 
k 

i=2 

Furthermore, the number of blocks Sn in a k-Stirling permutations of order n coincides with 
the number of left-right nodes in the corresponding {k + l)-ary increasing trees of order n, 

Sn = LRn. (17) 



Proof. Using the stated bijection we observe that a (j + 1) -child, I < j < k, corresponds to 
a j -ascent, since the step from the parent node v to the (j + 1) -child u corresponds to having 
recorded j times the label of the parent v and then another label w, with w > u. Similar 
considerations prove the results for j-descents and j-plateaux. The results for F„, 
then follow from dH)-® and ( fT6l) . 

Concerning the connection between blocks and left-right nodes we make the following 
observation. Starting with a (A; + l)-ary increasing tree, and inserting nodes one after another, 
we note that only a leftmost or rightmost child leads to a new block in the corresponding k- 
Stirling permutation. Hence, the number of left-right nodes is equal to the number of blocks, 
since we start with a single block 1^ and a single left-right node (the root). □ 
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Remark 3. Note that the number of leaves in {k + l)-ary increasing trees of order n > 2 
corresponds to the number of locally maximal substrings l'^ = I - ■ - l, i.e. substrings il'^j, 
with < i, j < I, for 2 < / < n, in /c-Stirling permutations of order n, which can also be 
seen from the bijection. 

Remark 4. In the case k = 2 we thus have the symmetric situation that X„ = I/„ i, F„ = 
L„ 3 and Z„ = L.„ 2, which by Theorem [8] below gives a new proof that X„, F„ and Z„ have 
the same distribution, and further are exchangeable, as shown by [6] and [18]. We see also 
that this will not hold for larger k, see for example Theorem|9l 

For k = 2, Theorem [T] gives a bijection between Stirling permutation and ternary increas- 
ing trees, while Janson [Ull gives a bijection with plane recursive trees of order n + 1. (These 
are related by a bijection given in Section HI) Next we will show that also for k > 2, there 
is a suitable family of generalized plane recursive trees that is closely related to /c-Stirling 
permutations. 

Definition 1. For k > 2, the family of k-plane recursive trees is specified by the degree- 
weight generating function Lp{t) = (l — {k — i.e. it is the family of generalized 
plane recursive trees with ip^ = 1, ci = k and C2 = —{k — 1). Explicitly, = ^ 11/=! ((^ ~ 
1)(/ — 1) + 1). Consequently, by (fT3l) and (fT4l) . the generating function T{z) and the total 
weight Tn are given by 



For k = 2, these are the plane recursive trees in Example [H 

Remark 5. We did not succeed in finding a bijective correspondence between /c-Stirling 
permutations and A;-plane recursive trees, in the case of /c > 2, generalizing the bijection in 
IfTSlI for k = 2, since for A; > 2 it seems difficult to obtain a combinatorial interpretation of 
the weights of the trees. We leave this as an open problem. However, the distribution of the 
leaves still coincides with the distribution of the number of ascents or descents. 

Theorem 3. The number ( total weight) of k-plane recursive trees of order n + 1 equals 
the number of k-Stirling permutations of order n, Tn+i = Qn- Moreover, the distribution 
of the number L^+i of leaves of k-plane recursive trees of order n + 1 coincides with the 
distribution of the number X„ of ascents (descents) of k-Stirling permutations of order n. 

Proof. The first part is already shown. 

The second part is trivial for n = 1, with one leaf and one ascent. We proceed by induction, 

and suppose that the relation is true for n: Ln+i = X„. We observe that adding the new node 
labelled n + 2 to a leaf does not change the number of leaves, whereas adding the new 




n 



k-l 
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node at any other place gives rise to a new leaf. Further, by the formula for p{v) above 

with d(v)=0, the probability of adding node n + 2 to a given leaf in a tree of order n + 1 is 

1 

p{v) = Tij^ — = ^^^V]-- Hence, conditioned on the number of leaves L^+i being m, we 

k-l 

have Ln+2 = m or m + 1 with 

P(L„+2 = m I L„+i = m) = — . (18) 

kn + 1 

Similarly, when adding a string (n + l)'^ to a fc-Stirling permutation of order n, we will 
always create a new ascent, and we will destroy one if and only if we add the string at an 
ascent. Since there are kn + 1 gaps where the new string can be added, conditioned on the 
number X„ of ascents being m, we have Xn+i = m or m + 1 with P(X„+2 = m \ Xn+i = 

m) = -j^^^- This is the same relation as (fTSl ), and thus Ln+2 = Xn+i, which verifies the 
induction step. □ 

Remark 6. The distribution of the number of leaves is fairly well studied. Let T(z, v) = 
J2n>i Tn,m^'v^ dcuotc the bivariate generating function of the number of A; -plane recursive 
trees having exactly m leaves, also encoding the number fc-Stirling permutations of order 
n — 1 having m descents. Bergeron et al. [3| determined the generating function T{z, v) by 
the implicit equation 

dt 



'o {v-l)^Q + ^(t) 

Note that the implicit equation is true for a much larger class of increasing trees; moreover 
one may derive the normal limit of the number of leaves from the implicit equation above, 
see 0. 

3.2. A;-bundled increasing trees and A;-bundled Stirling permutations. 

Definition 2. For A; > 0, the family of {k + l)-bundled increasing trees is specified by the 
degree-weight generating function Lp{t) = j^^z^jk+r, i.e. it is the family of generalized plane 

recursive trees with ipo = 1, ci = k + 2 and C2 = —1. Explicitly, (fj = {^^'')- Consequently, 
by (fT3l) and ([T4)) . the generating function T(z) and the total weight T„ are given by 

n-l 

T{z) = l-{l-{k + 2)^)^, T„ = l[{l{k + 2) - 1), 

1=1 

a = k+1, and the the probability attaching to node v at step i + 1 is given hy p{v) = '^^.^^^^^^ . 

Remark 7. One may think of {k + 1) -bundled increasing trees of order n as consisting 
of a root node labelled 1 which has k + 1 positions, with a (possibly empty) sequence of 
labelled (k + l)-bundled increasing trees attached to each position (with disjoint sets of 
labels, forming a partition of {2, . . . ,n}). Equivalently, one may think of each node as 
having k separation walls, which can be regarded as a special type of edges. 
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Note that the 1-bundled increasing trees are just ordinary plane recursive trees, cf. Exam- 
ple [3l and that the bijection stated below also holds for this case, which corresponds to the 
resultof Ea that;B„(l) = % = Q„_i (2), since obviously Q„(0) ^ Q„_i (2) by relabelling. 

Theorem 4. The family Bn = Bnik + 1) of{k + 1) -bundled increasing trees of order n is in 
a natural bijection with k-bundled Stirling permutations, Bn{k + 1) = Qn{k). 

Proof. We proceed as before using a depth-first walk. We label each auxilliary separation 
wall of a node labelled v by the label of the node v. Moreover, we label any (proper) edge by 
the label of the child. Hence, at any time, any node has at least k outgoing edges, thinking 
of the walls as a special type of edges. Now we perform the depth-first walk and code 
the fc-bundled increasing tree by the sequence of the labels visited on the edges, under the 
additional rule that a label on a separation wall only contributes once. Since every proper 
edge is traversed twice, and every label except 1 occurs on exactly one proper edge, a (A; + 1)- 
bundled increasing tree of order n is encoded by a string of ( A; + 2) (n — 1) + /c integers, where 
each of the labels 2, . . . , n appears exactly k + 2 times and label 1 appears k times. In other 
words, the code is a permutation of the multiset {l'^, 2^^"^, . . . , n''^'^}. Note that for each 
^, 1 < ^ < ^, the elements occurring between the two occurrences of i are larger than z, 
since we can only visit nodes with higher labels. Hence the code is a A;-bundled Stirling 
permutation. Moreover, adding a new node n + 1 at one of the {k + 2){n — \) + k + \ 
possible places corresponds to inserting the [k + 2)-tuple {n + in the code, at one of 
[k + 2){n — 1) + k + 1 possible places. This shows that the code determines the {k + 1)- 
bundled increasing tree uniquely and that the coding is a bijection. See Figure [2] for an 
illustration. □ 




Figure 2. The two 2-bundled increasing trees of order 2 encoded by 2221, 
1222; Three 2-bundled increasing trees of order 3 encoded by the sequences 
2333221, 3331222 and 3332221. 

Next we relate the distribution of ascents, descents and plateaux in /c -bundled Stirling 
permutations with certain parameters in (fc + 1) -bundled increasing trees. In order to do so 
we introduce three parameters for a {k -\- 1) -bundled increasing tree r. The parameter Ba 
counts the number of ascents in the bundles of r, plus the number of non-empty bundles, 
plus 1 if the first bundle of the root is empty, where an ascent in a bundle occurs if the root of 
a subtree is smaller then the root of the next subtree, going from left to right. The parameter 
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Be counts the number descents in the bundles of r, plus the number of non-empty bundles, 
plus 1 if the last bundle of the root is empty, where a descent in a bundle occurs if the root 
of subtree is larger then its neighbour. The number Be counts the number of empty bundles 
of the nodes with labels larger than one plus the number of empty inner bundles of the root. 
With these definitions, the following correspondences are straightforward. 

Theorem 5. Under the bijection in Theorem^ the numbers of ascents, descents and plateaux 
in a k-bundled Stirling permutation of order n coincide with the parameters Ba, Bo and Be 
in a {k + 1) -bundled increasing tree of order n. 

Remark 8. Note that the number of leaves in (A; + 1) -bundled increasing trees of order n 
corresponds to the number of sequences of the form l^^'^ = I ■ ■ - l, with 2 < Z < n, in 
A; -bundled Stirling permutations of order n, as for fc-ary increasing trees. Moreover, the 
parameter "number of descendants of node j" in a (A; + l)-bundled increasing tree of order 
n, with 2 < j < n, counts the number of different entries / with j < I < n between the first 
and the last occurrence of j in the corresponding /c -bundled Stirling permutation of order n. 

4. Further bijections 

The bijections of Theorem [T] (with k = 2) and niSH (or Theorem |4] with k = 0) imply a 
bijection between ordinary plane recursive trees of order n + 1 and ternary increasing trees 
of order n, using the connections to 2-Stirling permutations. In the following we will give 
two direct bijections, which both encompass this bijection between plane recursive trees and 
ternary increasing trees. 

First we give a bijection between sequences of fc-bundled increasing trees and {k + 2)-ary 
increasing trees, which for k = 1 just gives the desired bijection. 

Let SEQ(i3)„ = SEQ(i3)„(A;) denote the family of sequences of A;-bundled increasing 
trees with total order n, labelled with disjoint sets of labels forming a partition of {1, ... ,n}. 
(Note that our notation slightly abuses the common sequence notation SEQ of combinatorial 
objects, since we also assume properly distributed labels.) 

Remark 9. By introducing a new root labelled 0, connecting all roots of the sequence with 
the new root, and performing a proper relabelling, SEQ(i3)n is in bijection with the family 
of increasing plane trees of order n + 1 where each node except the root is A;-bundled as in 
Definition^ (Equivalently, SEQ(B)n is in bijection with the family of A;-bundled increasing 
trees of order n + 1 where the root has only the first bundle non-empty.) 

Theorem 6. The family SEQ(i3)„ = SEQ(i3)„(A;) of sequences of k-bundled increasing 
trees of total order n is in bijection with An{k + 2), the family of{k + 2)-ary increasing trees 
of order n: SEQ(-B)„(A;) = An{k + 2). 

Remark 10. Recall that 1 -bundled increasing trees are exactly plane recursive trees. More- 
over, in the case of A; = 1, the bijection in Remark |9] is the standard bijection between 
sequences of plane recursive trees of total order n and plane recursive trees of order n + 1; 
hence SEQ(i3)„(l) = See Figure[3]for an illustration. It is easily seen that, for k = 1, 
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the bijection Tn+i = SEQ(i3)„(l) = ^n(3) constructed in the proof below yields the corre- 
spondence between the two bijections ^„(3) = Q„(2) in Theorem [T] and T„+i = Q„(2) in 
[HI or-B„+i(l) = Qn+M = Q„(2) inTheoremia 

Proof. We use a recursive construction, see Figure [3] For a given sequence of A; -bundled 
increasing trees, we choose in the first step the tree of the sequence with node labelled 1: this 
node is going be the root of the (k + 2)-ary increasing tree. Since a (k + 2)-ary increasing 
trees has k+2 (possibly empty) subtrees Si, ... , going from left to right, we proceed as 
follows. The sequence of A; -bundled increasing trees to the left of the tree with root 1 forms 
(recursively) the subtree Si, conversely the sequence of A; -bundled increasing trees to the 
right of the tree with root 1 forms the subtree 5^+2. The k bundles, possibly empty, attached 
to the tree with root labelled 1, form the subtrees 5*2, ... , 5^+1 of the (k + 2)-ary increasing 
tree. Now we can proceed recursively, since the bundles are themselves just sequences of 
/c -bundled increasing trees. 

Conversely, starting with a (A; + 2)-ary increasing tree of order n, we recursively build 
a sequence of A; -bundled increasing trees as follows. In the first step we build a tree with 
root node labelled 1 . The sequence to the left of the tree with root labelled 1 is built from 
the subtree Si of the (A; + 2)-ary increasing tree of order n, the sequence on the right from 
the subtree Sk+2, and the k bundles are built from the subtrees 5*2, ... , 5*^+1. We proceed 
recursively until the sequence is constructed. Note that during this process, we connect any 
leftmost or rightmost child of a node v to the same parent as f . □ 




Figure 3. A sequence of 1 -bundled increasing trees of order 10, or equiv- 
alently a plane recursive tree of order 11, and the corresponding ternary in- 
creasing tree of order 10. 

Next we consider a bijection between A; -bundled increasing trees and so-called 
increasing trees. The family of fc+2-increasing trees consists of modified {k + 2)-ary 
increasing trees: any node except the root of a fc+2-increasing tree has A; + 2 labelled 
positions where children may be attached, whereas the root has only k positions (and thus 
outdegree bounded by A;). Note that for k = 1, the root has a single child and that chopping 
off the root yields a simple bijection between Fi 3 -increasing trees of order n + 1 and ternary 
increasing trees of order n. Thus the statement below implies for A: = 1 a bijection between 
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ternary increasing trees and plane recursive trees, An{S) = i3„+i(l) = which is just 

the bijection discussed in Remark [TOl 

Theorem 7. The family Tn = ^nik) of Fk^k+2-increasing trees of order n, is in bijection 
with the family of k-bundled increasing trees of order n, J-'n{k) = Bn{k), k > 1. 

Proof. For a given A; -bundled increasing tree of order n, we simply apply k times the bi- 
jection between sequences of /c-bundled increasing trees and (k + 2)-ary increasing trees to 
the k bundles attached to the root and the k positions of the root of the fc+2-increasing 
tree. □ 



Remark 11. To give an overview, we have provided the following bijections in Theo- 
rems [B II [6] and |71 



It is also possibly to give bijections Qn{k) = SEQ(i3)„(A; - 1) and Qn{k) = J^n{.k + 1), by 
simple modifications of the stated bijections. 

Remark 12. The families SEQ(i3)(A;) of sequences A; -bundled increasing trees and !F{k) of 
-Ffc,fe+2-increasing trees are non-standard in the sense that they are not part of the character- 
ization given by Panholzer and Prodinger \2T\. However, the counting problem concerning 
such tree families can be treated in a general manner, which will be discussed elsewhere. 

5. The distribution of j-ascents, j-descents and j-plateaux 

We are interested in the joint asymptotic distribution of j-ascents X„ j, j-descents Yn^j 
and j-plateaux Znj in a /c-Stirling permutations of order n, or equivalently in the joint dis- 
tribution of j-children Dnj and j-leaves Lnj+i in [k + l)-ary increasing trees of order n. 
Following Janson [I18II we use a (generalized) Polya urn model, see HISII . 

5.1. An urn model for the exterior leaves. Since we already know from (fT5l) that n — 
Dnj = Lnj, we can restrict ourselves to the study of the exterior nodes. We will use the 
following urn model. 

Urn I. Consider an urn with balls of A; + 1 colours, and let (L„ i, . . . , L„ fc+i) be the number 
of balls of each colour at time n. At each time step, draw one ball at random from the urn, 
discard it, and add one new ball of each colour. Start with {Li i, . . . , Li^^+i) = (1, 1, • • • , !)• 
Note that the vector (L„ i, . . . , Ln^k+i) exactly coincides (in distribution) with the numbers 
of the exterior nodes of types 1, . . . , /c + 1 in a random {k + l)-ary increasing tree, see 
Section [Xn 




Urn IJ is completely symmetric in the k + 1 colours, and we thus immediately see the 
following. 
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Theorem 8. For each n > \, the distribution of (Ln.i, • • • , -^n,fc+i) is exchangeable, i.e., 
invariant under any permutation of the k + 1 variables. 

It is customary and convenient to formulate generalized Polya urns using drawings with 
replacement. In the case of Urn HI we thus restate the description above and say instead that 
we draw a ball and replace it together with one ball each of the k other colours. In other 
words, Urn HI is described by the (/c + 1) x (/;: + 1) replacement matrix 

1 
1 



A — {1 — 5ij)i<ij<fc+i 



where 5ij denotes the Kronecker delta. 



1 

1 



1 
1 




1 

1 

M 



1 





1 
1 



1 



1 



1 
1 

0/ 



+ 1, have the same mean, 



5.2. Means. By Theorem [8l the variables L„j, j = 
and since their sum is A;n + 1 by (fT6l) . we see that they each have mean By (fTSi) and 

Theorem [2l we obtain the following exact formulas for the means. 

Theorem 9. The following hold, for n > 1 and k > 1: 

kn + 1 



EL 



ED. 



EX 



k+1 
n — 1 

k+V 
EY„ 



n — 1 



EZ 



k + 1 
kn + 1 



k+1 



l<j<k + l, 
l<J<k + l, 
l<J<k, 
l<3<k-l, 



EX„ = EK 



EZ, 



(k 



kn + 1 

k + 1' 
kn + 1 

^ k+1 ' 



5.3. Asymptotic distribution of j-ascents, j-descents and j-plateaux. We use the urn 

model Urn H] to obtain asymptotic normality. We begin with a general result. 

Theorem 10. Consider an urn with balls of q > 2 colours, where at each step one ball is 
drawn at random and discarded, and one ball of each colour is added. IfNnj is the number 
of balls of colour j after n steps, then, for any initial values Nq^i, . . . , N^^q, 

K --2.-1 



n 



V 
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jointly for j = 1, . . . , g, where Q are jointly normal random variables with means and 
(co)variances 

q - 1 

Note that Yli Cj = 0' for example because A^„j is deterministic. 

Proof. This urn has replacement matrix A = (1 — 6ij)1j^^ = (l)^j=i — Since the rank 
1 matrix {1)1 .^-^ has one eigenvalue q and q — 1 eigenvalues 0, A has largest eigenvalue 



•1 !•■ 



Xi = q — 1 and g — 1 eigenvalues —1. Theorem 3.22 in [1151 applies, with vi . . . , 

and shows joint convergence (in distribution) to normal variables Q with mean and a certain 
covariance matrix S. The formula for S in [JSi, Theorem 3.22] is complicated, so we use 
ifTSl Lemma 5.4], with a = {ai)l = {1)1 and m = Ai = g — 1, which yields S = mS/, with 
E, defined in [15, (2.15)]. Here = (^i.^O^Li = (1 - = {i^,jii,l)ll^^, vu = 1/q 

and, using the symmetry, B = {bij)jj^i with 6j j = ^ and bij = i ^ j; hence 
B = + Further, P/ is the projection onto the eigenspace of A for the eigenvalue 
—1, and thus P/ = lizli^zd Further, on this eigenspace A = —I and thus B = ^I, and ifTSl 
(2.15)] yields, noting that all involved matrices are symmetric and commute, 

= / Pie'^Be'^' P'je-^^' ds = -Pi / e^^-^-^"-^)^ = -Pi 

Jo q Jo + 

{q-l)I-A ( q6,,, - 1 



Recalling that E = mE/ = (g — 1)E/, we obtain the result. □ 

Remark 13. Similar calculations show, more generally, that if we at each step add a fixed 
number Si balls of colour % = 1, . . . , g, independently of the colour of the drawn and 

discarded ball, then n^^/^ (A^n.i — ^^^Siu) d, jointly, where d are jointly normal 
variables with means and 



EiSi + i\EiSi {Eisi)\ 

As an example, the numbers X„, F„ and Z„ of ascents, descents and plateaux in a random k- 
Stirling permutation can be seen as such an urn with replacement vector (1, 1, A; — 1), which 
yields an alternative proof of the limit distribution found above for them. 

We apply Theorem [TOl with q = k + 1, to Urn H] and obtain using (fTSl ) and Theorem [2] the 
following. 

Theorem 11. Let k > 1 and let Q, j = 1, . . . ,k + 1, be jointly normal random variables 
with means and (co)variances 

cov(o, 0) = + 1)*^^..- - 1)' 1 < ^, J < A: + 1, 
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fc2 



in particular Yai{Q) = (fc_|_i)2(fc_(_2) - ^^te that this implies Xlj=i = 0. Then, the following 
holds, jointly for all variables. 



r . 








\fn 




D 


~ fc+T'"' 












k+l"' 










Y 




(d) 
v^/ ^ 














































k(k-l) 
k+1^ 


{d) 


Jn 



0, 
-0, 

Vj ■= -Cj, 



l<j<k + l, 
l<j<k + l, 
l<J<k, 
l<J<k, 
l<J<k-l, 



V 



Ck+l, 
k 



i=2 



A simple calculation shows that the covariance matrix of i], Q is (cf. Remark [T3l) 



fc2 



(fc+l)2(fc+2) 



(fe+l)2(fc+2) 
fc2 



fc(A;-l) 



(fc+l)2(fc+2) 
k(k-l) 



\ 



(fc+l)2(fc+2) (fc+l)2(fc+2) (fc+l)2(fc+2) 
, k{k-l) k{k-l) 2fc(fc-l) , 

\ (fc+l)2(fc+2) (fc+l)2(fc+2) (fc+l)2(fc+2) / 

For k = 2, this yields the univariate limit theorems by Bona f6j| and the multivariate limit 
theorem by Janson BlBll for X„, Y^, Z„. 

For k = 1, the result for X„ or F„ reduces to the classical result on the asymptotics of the 
number of ascents or descents in a random permutation. (In this case Z„ = 0.) 

6. The distribution of the number of blocks 

The number Sn of blocks in a random fc-Stirling permutation is described by another urn 
model. 

Urn II. This urn has balls of two colours, black and white. At each time step, draw a ball at 
random from the urn, replace it and add k further balls: if the drawn ball was black, add k 
black balls; if the drawn ball was white, add 1 white ball and k — 1 black. Let Bn and Wn be 
the numbers of black and white balls in the urn at time n, and start with Wn = 2, _B„ = A; — 1. 
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We thus have Bn + Wn = kn + 1 balls in the urn at time n, and it is easily seen that the 
number of white balls can be interpreted as the number of gaps between the blocks, or first or 
last, in a random /c-Stirling permutation of order n, i.e. as the number of gaps where addition 
of a string (n + 1)'^ create a new block. This is one more than the number of blocks, and thus 
we have the equality in distribution 

Sn = W^- 1. (19) 

Urn In] is thus a 2 x 2 generalized Polya urn with ball replacement matrix M = {^'lii). 
This urn model is a special case of the triangular 2x2 urn models analysed in detail by Jan- 
son [|17 J . where the asymptotic distribution is given. The special case of balanced triangular 
2x2 urn models was also studied by Flajolet et al. |[T2l|. (An urn is called balanced if the 
total number of added balls is constant, independently of the observed color.) For the special 
case treated here we can add the exact distribution using the tree representation, the moments 
of Sn, and almost sure convergence. 

Theorem 12. The probability mass function of the random variable Sn counting the number 
of blocks in a random k-Stirling permutation of order n is given by 

The binomial moments E('^"^'^'^) are given by the explicit formula 



E 



The random variable Sn '■= ^^i+i i^n + I) is a positive martingale and converges almost 



'1 



surely to a limit (, i.e. Sn ^""^'\ further 



The limits ( and ( can be specified by the moments 

E(a = (r + l)!^^|^^, r>0. 

Further, ( has a density function f{x) that can be written as f{x) = r(i)x^^(7(x^^), x > 0, 
where g is the density function of a positive ^-stable distribution with Laplace transform 
e~'^^'^ ; it is thus given by the series expansion 
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Remark 14. The simple structure of the binomial moments and the almost sure convergence 
is actually true for all balanced triangular urns of the form M = ( ° ) , < a < /5, which 
is easily seen to be true by extending the martingale arguments to thie general case. 

Proof. We use three different approaches to study the block structure Sn in A;-Stirling permu- 
tations or equivalently the number of left-right edges Li?„ in {k + l)-ary increasing trees, see 
(fTTl ). In order to obtain the explicit results for the probability distribution of Sn, we analyze 
LRn- We can use the tree decomposition (flOl) in order to obtain the differential equation 

^T(z, v) = v{l + T{z, v)f{l + T{z)f-\ T(0, v) = 0, 

where T(z, v) = Yl,n>i Yl,m.>i ^{^n = ^}Tn^'v™' denotes the bivariate generating function 
of the numbers f'{Sn = rn}Tn, and T{z) = T{z, 1) is the generating function of total 
weights of {k + l)-ary increasing trees. By ExampleSl 1 + T{z) = (1 — kz)^^^^. Solving 
the differential equation and adapting to the initial condition gives the solution 

^ i-v{i-ii-kzy/>') ~ ^" 

Extraction of coefficients gives then the stated result for the probability mass function. More- 
over, the stated binomial moments may be obtained from the generating function by extract- 
ing coefficients. 



\ r J Tn I — w ^ I — w 



For the almost sure convergence we proceed as follows. Let Wn = S'„ + 1 be the number 
of gaps between blocks, or, equivalently, the number of white balls in Urn|nl see (fT9l ). Let JF„ 
denote the a-field generated by the first n steps. Moreover denote by A„ = Wn — Wn-i = 
Sn — Sn-1 G {0, 1} the increment at step n. We have 

E{Wn I J^n-l) = HWn-1 + A„ | Tn-l) = Wn-l + E(A„ | J^n-l) ■ 

Since the probability that a new white ball is generated at step n is proportional to the number 
of existing white balls (at step n — 1), we obtain further 

Hence, 

E{Sn I = [7;^|!:|^~||!y ._, = Sn-u n>2. 

f"-i+fe) k[n - 1) + 1 

V n— 1 / 

Hence, Sn is a martingale. Since it is a positive martingale, it converges almost surely to 
a limit (. By the well-known asymptotic formula ("^") ~ n"/r(a + 1), for any fixed real 

~ r (provided Sn oo), and thus Sn — — ^ C can also be written 

n-'/'Sn ^ C 
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More generally, we similarly have for any positive integer r 

'Sn + r 



E 



r 



r / \ r — 1 J k{n — 1) + 1 
+1 

k 



Sn-i + r\ n - 1 + ^ 



r / n - 1 + i 



I ' k 

Hence, ('^"^^) is a martingale, which also leads to the stated result for the moments 



fn-l+i 



ln-l + 

n-1 



in an alternative way, using the recurrence relation for the unconditional expectation. 
Letting n ^ oo in the moment formula yields 

fSn + r\^ ^ (r + l)4^i±i|-n^'A 

\ r ) ^ ^ (n-i+i) ^ ^r(l + r±i) 

which leads to ES*^ ~ {r+l)\r^^^^^^-rC^'^ . Hence all moments of n^^/^Sn converge, and the 

limits must be the moments of Q. Letting r ^ oo, we see that the moments do not grow too 
fast so that the moment generating function Ee*'' is finite for all t, and thus the distribution is 
determined by the moments. 

Finally, we use the general results for urn models in Janson ifTTl . First, [fTTl Theorem 

l-3(v)] yields the convergence Wn C iii distribution, and [fTTl Theorem 1.7] yields the 
moments of ^ that we just derived in a different way; note however that [T7, Theorem 1.7] 
yields the formula above also for non-integer r > 0, with the standard interpretation (r + 
1)! = r(r + 2). Furthermore, [,17. Theorem 1.8] shows that C has a density and gives the 
explicit formulas stated above. □ 

7. The sizes of the blocks 

Recall that every block in a fc-Stirling permutation begins and ends with the same label, 
which we can regard as a label of the block. We order the blocks in the block decomposition 
as /Ci, . . . , /Cs according to this label (where s is the number of blocks); thus /Ci is the block 
extending from the first 1 to the last, /C2 is the block formed by the smallest label not in /Ci, 
and so on. We also let Ki := |/Cj| denote the size of the ith block in this order, and put 
^ . = 0, ^. = for i > s. 

Alternatively, we may order the blocks according to decreasing size. We let Ki > K2 > 
. . . be the sizes of the blocks in this order, again with Ki = for i > s. Thus, {Ki)f is the 
decreasing rearrangement of (Ki)"^. 

For a random /^-Stirling permutation of order n, we use the notations /C„,j, Kn,i and Kn,i. 
Note that Y,i Kn,i = Zlj ^n,i = kn. 

To study these sizes we introduce another urn model. Consider first an urn with balls of 
two colours, Kn,i — 1 white balls representing the gaps inside the block /C„^i and nk+2—Kn^i 
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black balls representing the gaps outside. Adding the string (n+ 1)^ at one of the gaps inside 
lCn,i means increasing Kn,i by k, and adding it outside means keeping Kn,i unchanged; 
hence this is a Poly a urn of the original type considered by Eggenberger and Poly a [fTTTl . 
Il28l . where we draw a ball at random and replace it together with k balls of the same colour. 
We start with Ki i = k, and thus k — 1 white and 2 black balls. 

Next, let us study the second block, /C„,2- At the first n where this is non-empty, we have 
k + 2 gaps outside the first block /C„ i, /c — 1 of them in /C„ 2 and 3 outside both blocks. Let 
us now ignore the first block and consider an urn with Kn.2 — 1 white balls representing the 
gaps in /C„ 2 and black balls representing the gaps outside both /C„ 1 and /C„ 2- The balls in 
this urn are drawn at random times (when we do not add to a gap in /C„ 1), but when they are 
drawn, the urn behaves exactly as for ICn,i'- we replace the drawn ball together with k of the 
same colour. 

The same argument applies to /C„ „ for any m > 2; if we ignore the preceding blocks and 
additions to them, we have the same Poly a urn again, but now started with m + 1 black balls, 
representing the gaps outside the first m blocks. We hence make the following definition. 

Urn III. This is the standard Pdlya urn with balls of two colours and where each drawn ball 
is replaced together with k balls of the same colour. Let Urn lllll „ be the version where we 
start with k — 1 white and m + 1 black balls, and let WV,™ and Bj^ „i denote the numbers of 
white and black balls after — 1 draws, when the urn contains W^^jn + -Biv.m = kN + m 
balls. 

We can thus identify (with the urns Urn lllll . Um lllll . . . . independent), recalling that the 
balls in urn UmlIIlL+i correspond to the black balls in urn UrnllllL. 

Kn,l = Wn,l + 1, 
Kn,2 = Wn,,2 + 1, 

and so on. 

Theorem 13. There exists a sequence of independent beta distributed random variables 
(3m ~ Beta(^, ^) such that 

K„,2, . . . ) ^ (A, (1 - /5i)/52, (1 - - ...). (20) 



with kN2 + 2 = Bn,i, 
with kN^ + 3 = 1, 



Proof. The basic limit theorem for Poly a urns says that, as — 00, 

WNm {as.) ^ „ fk — lm+l' 



I3m ~ Beta^- 



kN ' \ k ' k 
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and thus (as.)^ ^ _ (This is already in Polya [|28l for convergence in distribution. 

See, for example, [fT9l or ffTTl. Section 1 1].) Consequently, 

Kn.l Wn,l + 1 [a.s.) 

^ 



kn kn 

Kn,2 Bn,l Wn2,2 + 1 {a.s.) 



kn kn kN2 + 2 ' 
and so on. □ 



Note that both sides of (|20l ) are elements of V, the space of non-negative sequences {pi)f 
with J2iPi = 1; ^ can be seen as the space of probability distributions on N. The conver- 
gence in the proof above was componentwise, i.e. in the product topology, but it is well- 
known (and easy to verify) that on V, this topology is equivalent to the £^ -topology with 
the metric d{{pi), (p-)) = \Pi ~ Pi\^ ^nd also to the usual weak topology of probability 
distributions; hence the theorem holds for any of these topologies. 

Let Vi = Pi 11}=! (1 ~ f^j) be the elements of the limit sequence in (|20|) . and let (Vi)f° 
denote the decreasing rearrangements of them. The distribution of this random element of V 
is denoted PD(|, |), see Pitman and Yor ||27l or Bertoin [4J. 

Taking the decreasing rearrangement is a continuous operation on V, and thus we imme- 
diately obtain from Theorem [13] the following. 

Theorem 14. 

. . . ) ^ (VI, ^2, . . . ) ~ PD (l, I) . (21) 

kn \k k/ 

Corollary 1. The largest block size has the limit 

Kn,l {a.s.) ~ 

— > Vi = max Vi. 

kn i>i 

Remark 15. These results can be compared with the classical result that the lengths of the 
cycles in a random permutation, arranged in decreasing order and divided by the size of the 
permutation, converge (in distribution) to PD(1) = PD(0, 1), see e.g. |l2l Sections 5.5 and 
5.7]. 

Remark 16. For /c = 2 we obtain in Theorem [14] the limit distribution PD(i, i) which 
arises in other contexts too: it is the distribution of the sequence of excursion lengths in a 
Brownian bridge |[26l . |[25l . [[T], [|27l (for a related characterization for k > 2 see [27]) and it 
is the asymptotic distribution of the sizes of the tree components in a random mapping, see 
[[291 and [[2. It is an interesting problem to see whether there are more direct relations with 
these objects. 
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Schrodinger Institute in Vienna, 2008. 



24 s. janson, m. kuba, and a. panholzer 

References 

[1] DJ. Aldous and J. Pitman, Brownian bridge asymptotics for random mappings. Random Structures Al- 
gorithms 5, no. 4, 487-512, 1994. 

[2] R. Arratia, A. D. Barbour and S. Tavare, Logarithmic Combinatorial Structures: a Probabilistic Ap- 
proach, European Math. Soc, Ziirich, 2003. 

[3] F. Bergeron, R Flajolet and B. Salvy, Varieties of increasing trees. Lecture Notes in Computer Science 
581, 24-48, Springer, Berlin, 1992. 

[4] J. Bertoin, Random fragmentation and coagulation processes. Cambridge Univ. Press, Cambridge, 2006. 

[5] B. BoUobas and O. M. Riordan, Mathematical results on scale-free random graphs, in Handbook of 
Graphs and Networks, 1-34, Wiley-VCH, Weinheim, 2003. 

[6] M. Bona, Real zeros and normal distribution for statistics on Stirling permutations defined by Gessel and 
Stanley. Preprint, 2007. Online availaible at: http://arxiv.org/pcif/0708. 3223 .pdf 

[7] F. Brenti, Unimodal, log-concave, and Polya frequency sequences in combinatorics, Memc»j>iAmer.Ma?/i. 
Soc. 81, no. 413, 1989. 

[8] F. Brenti, Hilbert polynomials in combinatorics, J. Algebraic Combinatorics 7, 127-156, 1998. 
[91 D. Y. C. Chan, B. D. Hughes, A. S. Leong and W. J. Reed, Stochastically evolving networks. Physical 
Review E 68, 066124, 2003. 
[10] R. Dobrow and R. Smythe, Poisson approximations for functionals of random trees. Random Structures 

Algorithms 9, 79-92, 1996. 

[11] F. Eggenberger and G. Polya, Uber die Statistik verketteter Vorgange. Zeitschrift Angew. Math. Mech. 3, 
279-289, 1923. 

[12] P. Flajolet, P. Dumas, and V. Puyhaubert, Some exactly solvable models of urn process theory. Discrete 
Mathematics and Computer Science, Proceedings of Fourth Colloquium on Mathematics and Computer 
Science, Ph. Chassaing Editor., vol. AG, 59-1 18, 2006. 

[13] P. Flajolet & R. Sedgewick, Analytic Combinatorics. Cambridge Univ. Press, 2008, to appear. 

[14] I. Gessel and R. P. Stanley, StirUng polynomials. / Combin. Theory Sen A 24, no. 1, 24-33, 1978. 

[15] S. Janson, Functional limit theorems for multitype branching processes and generiiUzed Polya urns. Sto- 
chastic Processes Appl. 110, 177-245,2004. 

[16] S. Janson, Asymptotic degree distribution in random recursive trees. Random Structures Algorithms 26, 
69-83, 2005. 

[17] S. Janson, Limit theorems for triangular urn schemes. Probability Theory and Related Fields 134, 417- 
452, 2006. 

[18] S. Janson, Plane recursive trees, Stirling permutations and an urn model. Preprint, 2008. Online availaible 

at: http://arxiv.org/pdf/0803 . 112 9vl 
[19] N. L. Johnson and S. Kotz, Urn models and their application, Wiley, New York, 1977. 
[20] M. Kuba and A. Panholzer, On the degree distribution of the nodes in increasing trees, J. Combin. Theory 

SerAlU, 597-618, 2007. 

[21] A. Panholzer and H. Prodinger, The level of nodes in increasing trees revisited. Random Structures and 
Algorithms, 31, 203-226, 2007. 

[22] S. K. Park, The r-multipermutations, / Combin. Theory Sen A 67, no. 1, 44-71, 1994. 

[23] S. K. Park, Inverse descents of r-multipermutations. Discrete Math. 132, no. 1-3, 215-229, 1994. 

[24] S. K. Park, P-partitions and g-Stirling numbers. J. Combin. Theory Sen A 68, no. 1, 33-52, 1994. 

[25] M. Perman, J. Pitman and M. Yor, Size-biased sampling of Poisson point processes and excursions, 
Probab. Theory Related Fields 92, no. 1, 21-39, 1992. 

[26] J. Pitman &. M. Yor, Arcsine laws and interval partitions derived from a stable subordinator, Proc. London 
Math. Soc. (3) 65, no. 2, 326-356, 1992. 

[27] J. Pitman and M. Yor, The two-parameter Poisson-Dirichlet distribution derived from a stable subordina- 
tor. Ann. Probab. 25, no. 2, 855-900, 1997. 



STIRLING PERMUTATIONS, INCREASING TREES AND URN MODELS 



25 



[28] G. Polya, Sur quelques points de la theorie des probabilites. Ann. Inst. Poincare 1, 1 17-161, 1931. 

[29] V. E. Stepanov, Limit distributions of certain characteristics of random mappings (Russian), Tear. Vero- 
jatnost. i Primenen. 14, 639-653, 1969; English transl. Theory Probab. Appl. 14, 612-626, 1969. 

[30] J. Vitter and P. Flajolet, Average case analysis of algorithms and data structures, in Handbook of Theoret- 
ical Computer Science, 431-524, Elsevier, Amsterdam, 1990. 

[31] H. Wilf, Generatingfunctionology, 2nd ed., Academic Press, Boston, 1994. 

SVANTE JANSON, DEPARTMENT OF MATHEMATICS, UPPSALA UNIVERSITY, PO BOX 480, SE-751 06 

Uppsala, Sweden 

E-mail address: svante . jansonSmath .uu.se 
URL: http: //www. math. uu. se/~svante/ 

Markus Kuba, Institut fur Diskrete Mathematik und Geometrie, Technische Univer- 
siTAT WiEN, Wiedner Hauptstr. 8-10/104, 1040 Wien, Austria 

E-mail address: kuba@dmg . tuwien . ac . at 

Alois Panholzer, Institut fur Diskrete Mathematik und Geometrie, Technische Uni- 
versitat Wien, Wiedner Hauptstr. 8-10/104, 1040 Wien, Austria 
E-mail address: Alois . Panholzer@tuwien . ac . at 



